Vorticity cutoff in nonlinear photonic crystals 
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Using group theory arguments, we demonstrate that, unlike in homogeneous media, no symmetric 
vortices of arbitrary order can be generated in two-dimensional (2D) nonlinear systems possessing 
a discrete-point symmetry. The only condition needed is that the non-linearity term exclusively 
depends on the modulus of the field. In the particular case of 2D periodic systems, such as nonlinear 
photonic crystals or Bose-Einstein condensates in periodic potentials, it is shown that the realization 
of discrete symmetry forbids the existence of symmetric vortex solutions with vorticity higher than 
two. 



Vortices are particular higher-order stationary solu- 
tions present in many different nonlinear systems, rang- 
ing from fluid dynamics to photonics. A vortex is charac- 
terized by a typical phase dislocation determined by an 
integer number, that we refer to as vorticity (also known 
as winding- number, "topological charge" or even spin). 
An optical vortex with a rotationally invariant amplitude 
in a nonlinear Kerr medium, experimentally observed in 
homogeneous self-defocussing media 0, can be under- 
stood as an eigenmode of the equivalent rotationally in- 
variant waveguide generated by itself 0- Thus, a vortex 
appear as an object carrying well-defined angular mo- 
mentum: (pi — e*'^/(r). In this case, angular momentum 
and vorticity are the same integer number; a consequence 
of the continuous 0(2)-symmetry of the operator defining 
the equivalent waveguide. However, in systems such as 
2D nonlinear photonic crystals or Bose-Einstein conden- 
sates in 2D periodic traps this 0(2)-symmetry is replaced 
by a discrete point-symmetry. Angular momentum is 
no longer well defined and thus the angular-momentum- 
vorticity equivalence is lost. Nevertheless, optical vor- 
tices have been predicted to exist in 2D periodic photonic 
crystals 0,0] and in photonic crystal fibers and exper- 
imentally observed in optically-induced photonic lattices 
1^. Although these solutions cannot longer have well- 
defined angular momentum, certainly all of them present 
neat phase dislocations that can be characterized by an 
integer vorticity value. In this paper, we will prove how 
to re- interpret vorticity in terms of the rotational proper- 
ties of vortex solutions without resorting to the angular- 
momentum concept. As a result, severe restrictions on 
vorticity values will be found using group-theory argu- 
ments. 

Let us consider the following general nonlinear equa- 
tion for stationary states: 



[io + iNL(|0|)] v) = -£(l){x, y). 



(1) 



where Lq is a linear field-independent self-adjoint oper- 
ator (normally dependent on gradients and functions of 
the transverse coordinates) and Lnl(|0|) is the nonlinear 
field-dependent piece of the full operator acting on the 
field (j). This equation is valid for all type of 2D systems 



in which the nonlinearity depends on the field through 
its modulus. Many different systems can be modeled us- 
ing an equation that can be written in the form given by 
Eq.(^. We are interested in systems that, besides being 
described by Eq.JQ), are invariant under some discrete- 
symmetry group G: [L,G\ = (L = Lq + ^nl)- This 
means that we assume that all linear and nonlinear co- 
efficients appearing in the operators defining Eq.(^ are 
invariant under the G group. Our goal is to study the im- 
plications that the realization of discrete symmetry have 
on the characterization of vortex solutions of Eq. (QJ . 

The key concept in our approach is the so-called group 
self-consistency condition. This condition establishes 
that if a system described by Eq.l^ is invariant under 
some discrete-symmetry group G then any of its solu- 
tions either belongs to one representation of the group G 
or to one of its subgroups G' [G' C G) 7]. Note that the 
identity group is always a subgroup of any group and, 
therefore, asymmetric solutions also satisfy the group 
self-consistency condition 0. In this letter, however, we 
will focus on symmetric solutions exclusively. 

The elements of point-symmetry groups in a plane are 
rotations through integral multiples of 27r/n about some 
axis (called an n-fold rotation axis), reflections on a mir- 
ror plane containing the axis of rotation and combina- 
tions of both. Groups containing a n-fold rotation axis 
constitute the C„ groups. When, in combination with 
the n-fold rotation axis, these groups have mirror planes, 
one generates the so-called C„„ groups. In Fig. Q we 
give two examples of structures exhibiting Cq^, and Cs^ 
point-symmetries. 

We will prove next how vorticity is affected by the 
finite order of the n-fold rotation axis defining the C„ 
(or Cnv) group. In order to do so, we need first to 
properly characterize the different representations of a 
Cn group. Since C„ groups are abelian, its representa- 
tions are one-dimensional and given by a single scalar 
complex number (the character of the representation) 
0. This scalar is nothing but a root of unity of or- 
der n and thus the representations of the C„ group 
are given by { 1 , e 
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Figure 1: Two examples of structures invariant under 27r/n 
rotations plus specular reflections on the x and y axis; (a) 
6-fold rotation axis (Ce^ group), and (b) 8- fold rotation axis 
{Csv group). 
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Figure 2: Roots of unity diagrams displaying the representa- 
tions of; (a) Ce and (b) C3. 



exp(27ri/rt). In Fig. |21 we present, as an example, the 
construction of the roots of unity for the Cg (even n) 
and C3 (odd n) groups. Each representation can be la- 
beled by the natural number I and, when present, by its 
sign. We denote it by Vi^s (Z g N, s = ±). No sign 
is needed for the identity representation 2?o (' — 0) nor 
for I7„/2 {I = ?T./2, even n). A state belonging to repre- 
sentations with I ^ 0, n/2 can be written as \l,s) with 
< / < n/2 (if n is even) or < ^ < (n - l)/2 (if n 
is odd). When we act with a group operator G (rep- 
resenting a discrete rotation of angle 27r/n) on a func- 
tion belonging to a representation it transforms as 
G(j)i = e'(/)f, where I = si {s = ±,Z G N). Clearly too, 
G(j)o = 00 and G(f>„/2 = ^"^^0n/2 (even n). If there is no 
other symmetry involved, [L,G] = implies that every 
one-dimensional representation is characterized by a dif- 
ferent L-eigenvalue: L(f)i = —£i(j)i. Representations are 
thus non-degenerated. 

We proceed now to explicitly construct functions be- 
longing to Z 7^ representations. Let us consider the 
complex coordinate vector u ^ x + iy — re*^. Integer 
powers of u have well-defined transformation properties 

under a zn/n rotation: u e u . Therefore, we 

can easily construct a function in the representation 
of Cn as 



00 being a function in the "Dq representation of C„. 
Clearly, Gtpi — eVr- 

The representations of C„„ (discrete rotations plus re- 
flections) are easily obtained from those of C„ groups 
The existence of the extra symmetries provided by 
mirror reflections yields to degeneracies for I ^ rep- 
resentations. High-order states are now doubly degen- 
erated; they form pairs of complex-conjugated functions 
((/);,(/);*) with the same L-eigenvalue: L(\)i ~ —£i(j)i. Re- 
markable exceptions are the and X'„/2 representa- 
tions. Because of their different behavior under mirror 
reflections there are two distinct non-degenerated one- 
dimensional I — representations: |0; -I— f ) and |0; ). 

They transform differently under reflections with re- 
spect to X and y axis: i^^^^y |0; -I— 1-) = -I- [0; -I— 1-) and 



|0; 



The |0; state has maximal 



symmetry. Fundamental solitons belongs to this identity 
representation of Cnv In the same way, there are also 
two different non-degenerated one-dimensional represen- 
tations with I = n/2 (even n): \n/2; H — ) and |n/2; — h). 
The distinction is made by i^a^^y-reflections. In Fig. |3| 
we show the lowest order eigenfunctions of the spec- 
trum of a Cgti-invariant operator self-consistently gener- 
ated by a fundamental soliton solution (0fund = '?!'o,-i-+)- 
L = Lq + iNL(|0fund|) (see the flnal section of the pa- 
per for details on the physical system associated to L). 
We easily recognize, from lower to higher values of £, the 
0; -I— f ) self-consistent state (i.e., the fundamental soli- 
ton), the doubly-degenerated |1; ±) and |2; ±) states and 
the non-degenerated |3; H — ) state. The rest of the spec- 
trum, including continuum de-localized states, systemat- 
ically falls into the representations described above. 

Vorticity v can be defined as the integer variation (in 
27r units) that the phase of a complex field experiments 
under a 27r rotation around a rotation axis. Solutions 
with non-zero vorticity are called vortices of order v. 
They are characterized by their rotation axis, whose in- 
tersection with the 2D plane defines the vortex center, 
where their amplitude vanishes. If ^(r, 9) represents 
the phase of a complex vortex field of order v given by 
= |/„|e**, then $(r,6' -h 27r) - $(r,6') = 2ttv, where 
the polar coordinates are referred to a reference frame 
centered on the rotation axis. For systems enjoying a 2D 
point symmetry, this axis is naturally given by the n-fold 
rotation axis of the corresponding C„ (or Cnv) group. 

According to the group self-consistency condition, all 
symmetric solutions of Eq. in a system with C„ sym- 
metry have to lie on the representations of C„ or of any 
of its subgroups. Let us consider now a solution (jyj in 
the Vi^s representation of C„ given by Eq.(|2Jl. Its phase 

will be given by aig(j)i{r,9) = 16 + arg 0q ' (r, 0). Since 
0Q''(r, 9) is invariant under rotations, arg(/)|-(r, 9 + 2tt) = 
a,i-g<j)i{r,9) + 2-kI. Therefore we find the important rela- 
tion between the index representation and vorticity: 



#(u) = mVo^(u), 



(2) 



v = l. 



(3) 
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Figure 3: Lowest order eigenfunctions of a nonlinear operator 
L generated by a soliton solution in the identity (fundamental) 
representation of Ceu- The symmetry of the full operator is 
C(yv'- [L, Cei)] = 0. 

Vortices are thus solutions belonging to 2?/,s representa- 
tions with / ^ 0. There is, however, no vortex associ- 
ated to I = n/2 (even n). It can be proved that 4>n/2 
is a real field, so that its argument is a function that 
can only take the values or tt. More explicitly, from 

Eq.@, cos (^n6'/2 + arg(/)[|"/^^(r,6l)^, which has 

the phase behavior of alternating signs typical of a nodal 
soliton and not of a vortex 0- In Cnv, the behavior of 

the (t>n/2,-\ and 4>n/2. h functions is also of the nodal- 

soliton type, as one can check by observing the phase of 
the \3;+-) state in Fig. |3 

Let us summarize now our main conclusions. Firstly, if 
a system is invariant under a C„ or Cnv point-symmetry 
group, the solutions of Eq.lQJ belong to representations 
of these groups or of their corresponding subgroups. Sec- 
ondly, symmetric solutions of Eq.(^ are characterized 
by the representation index I, which has an upper bound 
fixed by the order of the group: I < n/2 (even n) and 
I < (n — l)/2 (odd n). Thirdly, the vorticity v of the 
vortex solutions of such a system has a cutoff due to Eq. 
(O and the upper bound for I: 

\v\<n/2 (even n) and |w|<(n-l)/2 (odd n). (4) 



Note that the group of continuous rotations on a plane 
can be understood as the limiting case 0{2) = lim„.^oo C„ 
and, thus, Eq.Q) correctly establishes the absence of a 
cutoff for it {\v\ < oo). 

When we deal with 2D periodic systems, the realiza- 
tion of discrete symmetry has particular features. This is 
a well-known problem in crystallography 0. The crys- 
tal structure is constructed according to a pattern that 
repeats itself to "tessellate" the 2D plane. Patterns, un- 
like objects, are invariant under translations (defined by 
the periodicity of the crystal) . This translation property, 
inherent to periodicity, determines that only certain col- 
lections of symmetry elements are possible for patterns. 
In other words, only patterns that exhibit a selected set of 
symmetries can "tessellate" the 2D plane. The important 
result for us here is that pattern periodicity establishes 
a restriction on the order of discrete rotations allowed in 
plane groups. Only n-fold rotations of order 2,3,4 and 6 
are permitted in a 2D periodic crystal 0- 

The previous group analysis has important implica- 
tions for 2D nonlinear periodic systems. According to 
the group self-consistent condition, if a periodic system 
described by Eq.(^ is invariant under a discrete group 
G, its stationary solutions cannot belong to representa- 
tions of groups with higher symmetry than G. If G' is 
the symmetry group of the solution, then G' C G. On 
the other hand, as seen above, the maximum rt-fold rota- 
tion symmetry compatible with periodicity is a sixth-fold 
rotation, which means that the maximum value for the 
order n of C„ and Cnv point-symmetry groups in 2D peri- 
odic systems is n = 6. Consequently, the point-symmetry 
group of a solution cannot exceed this order: n < 6. Since 
vorticity is restricted by the order of the point-symmetry 
group according to Eq. Q), we come up to the conclu- 
sion that in 2D nonlinear periodic systems of the type 
described by Eq.l^l vorticity has a strict bound: |w| < 2. 
Putting this into words, there are no vortices of order 
higher than two in 2D nonlinear periodic systems de- 
scribed by Eq.lpJ. 

In order to illustrate our previous theoretical results, 
we have numerically studied a realistic system, namely, 
a photonic crystal fiber (PCF). A PCF is a type of 2D 
photonic crystal consisting on a regular lattice of holes 
in silica (characterized by the hole radius a and the lat- 
tice period -or pitch- A) extending along the entire fiber 
length. When one considers that the silica response is 
nonlinear (nonlinearity represented by the nonlinear coef- 
ficient 7, defined in Ref.[S|), a PCF becomes a 2D nonlin- 
ear photonic crystal. The nonlinear propagation modes 
of a PCF for monochromatic illumination in the scalar 
approximation verify Eq.(^ with £ = ~0^, (3 being the 
mode propagation constant (see Among possible 

hole-distribution geometries we choose that based on a 
triangular lattice with Ce^ symmetry (see Fig. E^a)). 
The reason of our symmetry choice is simple. As proved 
before, the Cgt, group provides the highest vorticity solu- 
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Figure 4: Higher-order solitons for a periodic Cav PCF: (a)- 
(b) First- and second-order vortex pairs, |1; ±) and |2; ±); (c) 
nodal soliton of order three, |3; -I-—). 
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Figure 5: Same as in Fig. (|1J but in a PCF with defect. 

tions since it corresponds to the maximal point-symmetry 
achievable in a 2D nonlinear photonic crystal. In Fig. 
01 we find the first three (from lowest to highest value 
of /3^) higher-order solitons of a perfectly periodic PCF 
(without defect) calculated for the values a = 5 /im, 
A = 26 fim, and A = 1064 nm at 7 = 0.01. In Fig. El we 
present the same first three higher-order solitons but for 
a PCF with periodicity broken by the presence of a defect 
(absence of a hole). Note that in both cases the symme- 
try group is Cgt, and that, in agreement with our previous 
result, the maximum vorticity allowed is two. The soliton 



solution with Z = 3 is not a vortex. As predicted by group 
theory, it presents a binary phase structure (correspond- 
ing to a 1 3; -I — ) state) of the nodal-soliton type It 
is interesting to check the generality and accuracy of the 
group-theory approach using these numerical examples. 
The spectrum of higher-order soliton solutions is per- 
fectly explained by our previous group-theory arguments 
nevertheless the periodic (Fig. ^ and non-periodic (Fig. 
O photonic crystal structures present notable differences. 
Despite they share the same Cgu symmetry, a description 
in terms of weakly interacting localized fundamental soli- 
tons on lattice sites (the equivalent of the tight-binding 
approximation in solid state physics) 0| can only be valid 
in the perfectly periodic case. As it is apparent in Fig. 
121 this localization feature is clear in the amplitude and 
phase of vortex and nodal soliton solutions in the periodic 
PCF. However, single fundamental solitons are no longer 
recognizable in the vortex and nodal solitons of Fig. Odue 
to the presence of the periodicity-breaking defect. One 
can think of a situation of strongly interacting solitons 
causing the "tight-binding approximation" to stop being 
valid. Despite this fact, our main results concerning the 
nature of solutions and, more specifically, the restrictions 
on vorticity remain valid with complete generality. 
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